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Abstract 

Motivated by the recent measurements, we investigate B — > irir, Kir decay modes in 
the framework of QCD improved factorization, which was recently proposed by Beneke 
etal. We find that all the measured branching ratios are well accommodated in the 
reasonable parameter space except for B — ► K°7r°. We also discuss in detail the strong 
penguin contributions and the 0(a s ) corrections to the chirally enhanced terms. We find 
that the weak phase 7 lies in the region 120° < 7 < 240°, which is mainly constrained by 

B — > 7T~7r + . 

PACS Numbers: 13.25Hw,12.15Hh, 12.38Bx 



1 Introduction 

It is well known that the theoretical description of nonleptonic B decays is an extremely out- 
standing challenge, due to the nonperturbative nature of both initial and final mesons. A good 
understanding of the B nonleptonic decays, or at least a reliable estimation, is the prerequisite 
for extracting meaningful implications from experimental data and for testing the SM. In past 
years, some achievements have been performed toward the goal, for example, in Ref. [|T], 0, [|. 
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Recently, Beneke, Buchalla, Neubert and Sachrajda [|J have presented a promising factor- 
ization formula for the charmless nonleptonic B decays. The basic object in the calculation 
of B charmless nonleptonic decays is the hadronic matrix element (M 1 (p 1 )M 2 (p 2 )|(9j|i?(p)), 
where 0{ is the effective operator inducing the decay, M\ is the final meson absorbing the light 
spectator quark from B meson, and M 2 is another light meson flying fast from the b quark 
decay point as implied by Oi. The light spectator quark is translated softly to M\ and this 
effect would be taken to the nonperturbative form factor Ff 2 ~* Ml unless it undergoes a hard 
interaction. The quark pair, forming M 2 , ejected from b decay point carrying large energy 
of order of will involve hard interaction, since soft gluon with momentum of order Aqcd 
will decouple from the quark pair at leading order in AQCD/ m b i n the heavy quark limit. The 
essence of the argument of ||] can be summarized by the improved factorization formula 

(M 1 (p 1 )M 2 (p 2 )\O l \B(p)} = F B ^ Kh {Ml) f 1 dxT*{x)(j)M 2 {x) 

Jo 

+ / dxdydzT l II (x,y,z)(p Ml (x)(f)M2(y)4>B(,z), (1) 
Jo 

where <j)p(x) are the P meson's light-cone distribution amplitudes(DA). The hard amplitudes 
T^' 11 can be perturbatively expanded in a s (mj) and can be obtained from the calculations of the 
diagrams in Fig.l. It is interesting to note that T/ would be unity and T/ 7 would be absent at 
zeroth order of a s in the formula of Eq. (|TJ) , then the Naive Factorization would be reproduced. 
Another consequence of Eq.([l]) is that the final state interactions may be computable and 
appear to be the imaginary part of the hard scattering amplitudes. 

In this work, we extend the formalism to E> — > Kir decays and recalculate B — > ixn decays 
with electroweak penguin contributions. We also present detailed discussions about the strong 
penguin contributions and therefore we obtain the corrections to the chiral enhanced terms, 
which are found free of infrared divergence. We point out that there is large cancellation between 
the strong penguin hard scattering amplitudes and its contributions are small. Prospects of 
observing CP violation in those decay modes are also discussed. 
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2 Calculations 



First we begin with the weak effective Hamiltonian H e ff for the AB = 1 transitions as J 

Ueff ~ 71 

For convenience, we list below the operators in 7i e // for b — > q: 



/ 2 10 \ / 2 10 \ ■ 

VubV* q £ CiOf + y: a Oi + c g o g + v cb v; q £ c t o\ + £ ao, + c g o g . 

\i=l i=3 I \i=l i=3 ) . 



(2) 



Of = q a YLu a ■ u^^Lbp , 0% = q^Lup ■ u^^Lb a , 

0\ = q a YLc a ■ Cpj^Lbp , O c 2 = q a ^Lc p ■ Cp'j^Lba , 

3 = qal^Lba ■ Eg' q'pl^Lq'p , 4 = q a ^Lbp ■ £ g , q'^^Lq' a , 

5 = q a YLb a ■ J2 q ' q'plnRq'p , O e = q^Lbp ■ £ 3 , q'^^Rq' a , (3) 

C 7 = lq a YLb a ■ Eg' eqtq'pJuRq'p , 8 = ^q^Lbp ■ Eg' dq'q'pl^Rq'a > 

9 = ~q a ^Lb a ■ Eg' e^fipLip , Oi = ^q^Lbp ■ Eg' e q >q'^^Lq' a , 

9 = {g./^ 2 )m b d a (T^R{Xip/2)bpG^. 

Here q = d,s and (g'e{u, <i, s, c, 6}). a and /? are the ££7 (3) color indices and Afg, A = 1, 8 
are the Gell-Mann matrices; L and -R are the left- and right-handed projection operators with 
L — (1 — 75), i? = (1 + 75), and G^, denotes the gluonic field strength tensor. The Wilson 
coefficients evaluated at /1 = nib scale are0 

Ci = 1.082, C 2 = -0.185, 

C-i = 0.014, C 4 = -0.035, 

C 5 = 0.009, C e = -0.041, 

(4) 

C 7 = -0.002/137, C 8 = 0.054/137, 
C 9 = -1.292/137, C w = -0.262/137, 
C g = -0.143. 

After direct calculations, we get the hard scattering for the decay modes listed as follows 

Tp = % E P=u ,c V; q V pb Kiq^Lu) ® (u^Lb) + of (u lfl Lu) <g> (f^Lfe) 

+af (^W) ® (^6) + <{qiM) ® (q'^Lb) + a p 5 (q%Rq') ® (^L6) 
+a£(-2)(g J Rg') (g> (g'Lfc) + a?fe 9 ,(g' 7/iJ Rg') <g> (g 7 "L&) + (-2)(af + a 8a )(qRq') <g> (g'Lfc) 
+af|e (? '(g'7 M Lg / ) ® (g^Lfo) + (a? |e ? ' + a p 10a ){q^Lq') <g> (g^Lfo)], 

(5) 
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Figure 1: Order a s corrections to the hard scattering kernels T/ (a, b, c, d, e, f) and T/ 7 (g,h) 

where the symbol (g) denotes {M 1 M 2 \ j 2 ®ji\B) = (M 2 |j 2 |0}(Mi|j 1 |S}. The effective af's which 
contain next-to-leadingorder(NLO) coefficients and 0(a s ) hard scattering corrections are found 
to be 



'1,2 



°4 



p 



0, < = i = 3,5,7,8,9,10,8a, 10a. 
r C 2 u s C F 

C o C 



+ (C*4 + C^) Gm 2 (s/) + C g GM 2 ,g 
f=u 



C 6 a, C F 
^ C 5 a a C F 



12), 



CiG" M2 ( Sp ) + C 3 (G" A/2 ( Sg ) + G" M2 (s 6 )) + (C 4 + C 6 ) £ G' M2 ( S/ ) + C 9 



C 8 + ^, 



AT 



b 3 3 
(C 8 + C 10 ) £ -e/G' M2 ( S/ ) + C 9 -(e q G' M2 (s q ) + e b G' Ah (s b )) 

f=u 



N 4vr N - — 
ho - Cl0 + ^ + ^^ G9 ^, 



a 



a s Cp 



10a 



4tt N 



3 b 3 
(C 8 + C10)- e f G M 2 (s f ) + C 9 -(e q G M2 (s q ) + e b G M2 {s b )) 

f=u 



where q = d,s. q' = u, d, s and / = u,d, s, c, b. Cf = {N 2 — 1)/ (2 AT) and N = 3 is the number 
of colors. The internal quark mass in the penguin diagrams enters as Sf = m 2 Jm 2 . x = 1 — x 
and u — l — u. 

F = -vzfoJL-is + f^ + fn (7) 
nib 

i-i 1 — 2x 

fli 2 = J dxg(x)(j) M2 (x)i £((x)=3y— — lnx-3i7r, 

fS. = k jl±,„ t t t *^ M . (8) 



N /f 1 (0)Af| Jo z Jo x Jo * y 

r 1 2 

G M2 , g = ~ / dx-4>M 2 {x), (9) 
Jo X 

2 4 a r 1 r 1 

GM 2 (sq) = In h 4 / dx<pM 2 {x) \ du uu In [s„ — uux — ie] , (10) 

3 3 m?, Jo Jo 

g' M2 , 9 = -!lte\fi*M = -% (ii) 

G'M 2 {s q ) = - — In— + 3 / dx<p° M (x) f du uu\n[s q — uux — ie] , (12) 
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where <f)(x) and 0°(x) are the meson's leading- twist DA and twist-3 DA respectively. It should 
be noted that we have included 0(a s ) corrections to a,Q in Eq.@. Although the term in 
Eq.(||) is formally 1/M& suppressed, it is chirally enhanced by jj,p = Mp/(m q + m q i) and known 
to be important to interpret the CELO0 measurement. As a result the 0(a s ) correction to a§ 
would be the most important one among the corrections to Oj. We see that there are logarithm 
terms Infi/m^ appearing in Eqs.(|7|^[T2|), which is the result of one loop integration. If the 
scale fi is chosen small, the logarithm would be large and has to be resummed by using the 
renormalization group method. In this paper we choose \x = nib, then the logarithm disappeared 
and the resummation is not necessary. As a result, the effective coefficients af's are obtained 
to the order of a s (m&) corrections ( see also in Ref. 0). 

We realize that the contribution of the strong penguins depicted in fig.l.(e) and (f) to a Q 
could be reliably estimated without IR divergence. As an example, we show the contribution 
of Fig. 1(f) in the following. With the assignment of the vertex 8 a pifM 2 ^M 2 lb4 >0 { x ) /^N c to M 2 



and its constituents, we can get the hard amplitudes of Fig.l.(f) as 

atsC F f 1 o 3(1 -x)ml_ 



H f ~ ifu^hh-^-^- J Q dx<p ( x ) p 9i7/*(l-75)6i 



ft7 M (l-75)6i / dar0°(s). (13) 



o 

We can see that the end point IR divergence in l/k 2 (k 2 = (1 — x)m 2 ) is canceled by the term 
(1 — x) in the numerator and the amplitude is finite. For the amplitude of Fig.l.e, it is easy to 
note that the denominator k 2 of the gluon propagator is canceled by the quark loop and the 
integration of Jq 1 dxG(sf) is also finite itself. However, if all the external quarks are treated as 
free qurks at first, IR divergence will appear. In the case of free quarks, one can get the hard 
amplitudes of Fig.l.(f) as 

Tn — 



k 2 
ml 

k 2 



dily.il - 7 5 )&^7 M (1 - 75)% + da^l - 75)6^7^(1 + 75)% . (14) 



At this stage the quark pair qd is in color-singlet configuration. After Fierz rearrangement, 
one gets 



o 

771 

'dij^l - 7 5 )gi <8> gj7 M (l - l5)bj - 2dj(l + 7 5 )g, ® q~j(l - 75)^. (15) 



A- 2 

From the above equation we can see that Fig. 1(f) contributes to 04 and a 6 equally and its 
contribution is IR divergent when k 2 — > in free quark approach. Phenomenologically, one 
may have to treat k 2 as a parameter. In the framework employed here, the virtuality of the 
gluon is convoluted with the meson's DA. Furthermore, The NLO strong penguin contributions 
to 04 and terms are different. 

Finally, the chirally enhanced contribution of Fig. 1(g) and (h) to a% is cancelled between 
them. One can easily see this cancellation by putting both the leading-twist DA and twist-3 
DA 4>(x) and 4>°(x) to Fig. 1(g) and (h) and calculate these two diagrams. Because 4>°(x) gives 
the chirally enhanced contributions, one can easily see that these contributions are cancelled. 

With Eqs.([|)and (^|), we can write down the amplitudes of B — > txti and Kn decays 

+R n -(a>Q + a§ + a 8a )] - [a c 4 + a c w + a c Wa + R n (a c 6 + a c 8 + a 8a )]|, (16) 



M(B° d 



71%°) 



M{B- 



vrV 



M(B° d -±K°ir Q 



F . 



tU(M 2 B -Mt)F B ^(0)\\V cb \ x 



G 



3 3 1 
— a 2 + a 4 + -a 7 — —a 9 — t^io + a ioa 



i 

-a 



10 ~ "10a 



^lUMl - M*)F B ^(0)\\V cb \\ R b e~^ [< + a£ 



3 

+ -l-a 



- 2 [~ a 7 + ^° a 8 + a 9 + a lo)] 



R' b e-^ 



G, 



«2 - 2 K - <) 

r2\ 77B 



4 7 "9, 



i/*(M£ - M^)F^(0)(1 - A 2 )|K b | x 



(17) 



7T 



R' b e~^ 



Rk I a " — ~ a « 



- ( ag — -a 8 + a 8a 

r 2 \ jpB^K 



ia ) + -Oio 



o a 10 a 10a 



l 10a 



+ 

^jfUMl - M K )F^«(0)(1 - A 2 )|V C0 | x 
i^e"* 7 [< + < + R K (ag + ag + a 8a ) + a^, 



10 ~ "lOaJ 



[04 + ifo (ag + ag) + a 



10 1 "lOaJ 



G, 



t'/^M^ - M^)F^(0)(1 - A 2 )|V C0 | x 



(19) 



(20) 



M(B~ -> Ji 7T' 



# 6 e" 47 [a? + c# + ifo (a« + ag + a 8a ) + a? + a? J 
[a% + R K (ag + a 8 + a 8a ) + a^ + aj 0a ] | 
+ ^iU{Ml - M K )F B - K (0)(1 - ^)\V cb \ x 



R' h e~ ij 



<4 + 2 («9 " a?J 



+ 2 ( fl 9 ~ 



^M^l - M^)F B -(0)(1 - ^)\V cb \ x 



04 + Rk I «g — — ag + a 8a ) — o a io + a 



noa 



+ 



a\ + ag — -a 8 + a 8a 



-a 10 + a 10a 



(21) 



(22) 



Where R b 



l-A 2 /2 I V„„ 



A I V^t I 
l-A 2 /2 I V ci) I 



and 

the phase of V* b . A = = 0.2196. R P = 2/i P . 



V cb , V u d and V us are chosen to be real and 7 is 
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3 Numerical calculations and discussions of results 



In the numerical calculations we use || 

U = 0.133GW, f K = 0.158GeV, f B = 0.180GeV, 

r(B + ) = 1.65 x l(T 12 s, r(B°) = 1.56 x l(T 12 s, 

M B = 5.2792GW, M b = A.8GeV, M c = lAGeV, 

m u = 4.0MeV, M d = 9.0MeV, M s = 80MeV. 

For the leading- twist DA 4>(x) and the twist-3 DA <j)°(x) of K and ir, we use the well known 
asymptotic form of these DA[^, |T^| 

</> n>K (x) = 6x(l-x), < Ar (x) = l. (23) 



For B meson, the wave function is chosen as that used in |JTT|, [12 

4>b(x) = N B x 2 (l — x) 2 exp 



M 2 B x 2 
2u% 



(24) 



with ub = 0.4 GeV, and Nb is the normalization constant to make f$ dx(f>B(x) = 1. Here 
the decay constant in the wave function has been factored out. So the wave function can be 
normalized to 1. It is also necessary to note that (j) B {%) is strongly peaked around x = 0.1. 
This character is consistent with the observation of Heavy Quark Effective Theory that the 
wave function should be peaked around Aqcd/M B - With such choice, we find 

/ 1 rfx^ = 11 . 15; (25) 

JO X 

which is near to the argument |4[] in which Jq 1 dx(pB(x)/x = Mb/Xb = 17.56 with A# = 0.3GeV. 
We have used the unitarity of the CKM matrix V* q V ub + V* q V cb + V tq V t b = to decompose the 
amplitudes into terms containing V* q V U {, and V* q V c b, and 

\V ud \ = 1 - A 2 /2, \V ub /V cb \ = 0.085 ± 0.02 

\Vcb\ = 0.0395 ± 0.0017 \V US \ = A = 0.2196. 
We leave the CKM angle 7 as a free parameter. For the form factors, we use ^^(0) = 0.3 
and F B ^ K (0) = 1.13F B ^(0). 
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Numerical values for a\ (nn) and a%(nK) are presented in Table 1. It should be noted that 
ai(Kir) are generally different to ai(nir) and also change from case to case due to fU in the 
formulas of aj, where M 2 could be K or n. However, with our choice of parameters 

h. ~ iH (27) 

and the same DAs <f>K n {%), the a^Kir) ~ ai(mr). From Table. 1, we can find that all a\ 
develop strong phases due to hard strong scattering. Our a 2 is very different from that of 
|13|) [14| in both real and imaginary part because of the contribution of Fig.l.(g), and (h). So, 
theoretical predictions for the decays dominated by a 2 may be very different between Naive 
Factorization approach and QCD improved factorization approach. Numerically, we find that 
the 0(a s ) strong penguin contributions which collected in a^ and a§ are small because of the 
large cancellation between Fig.l.e and Fig.l.f. In detail, the strong penguin contributions to 
a 4 and a 6 are 



p a s Cp 



b 

C\Gm 2 ( s p) + C^{GM 2 { s q) + CrM 2 ( s b)) + (C*4 + Cq) ^ G M 2 ( s /) + C ' g G ' M 2 ,g 

f=u 



a s C F (-0.780 - 1.744*) + (0.858), p = u, 

x < (28) 

Aix N (-1.473 - 0.529?) + (0.858), p = c. 



p 

®6pen 



a s C F 
4>lV 



CiG' M2 (s p ) + C 3 (G'M 2 (s q ) + G' Ah (s b )) + (C 4 + C 6 ) £ GV 2 (s/) + C fl G" M2)fi 



a s C F (-0.780 - 1.299i) + (0.2145), p = u, 

— x \ (29) 

An N (-1.095 - 0.510i) + (0.2145), p = c. 

where the numbers in the brakets are the contibutions of Fig.l.e and Fig.l.f respectively. The 

cancellation in a$ is weaker than that in a^, since the contribution of Fig.l.f to a§ is small. The 

other diagrams will dominate the 0(a s ) hard scattering amplitudes. 

Now it is time to discuss branching ratios and CP asymmetries of B — > Kn and B — > nn 

in the QCD improved factorization approach. The branching ratio is given by 

Br{B -> Kn,nn) = t b / (lQnm B )\M(B -> Kn, nn)\ 2 s, (30) 

where s = 1/2 for B — > n°7r° mode, and s — 1 for the other decay modes. For the charged B 
meson decays, the direct CP asymmetry parameter is defined as 

A Mr = \M(B+ -> /)| 2 - \M{B- -> /)| 2 
CP - /)| 2 + -/)| 2 " 1 J 
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Table 1: The QCD coefficients af at NLO for renormalization scale \l = m;, ( in units of 10 4 
for ci3,...,aio). Results from different references are shown for comparison. 



Go 
3 



ours 
1.042+0.014i 
0.046-0.082i 
65.2+26.8i 



1.038+0.018i 
0.082-0.080i 
40+20i 



m 

1.05 
0.053 
18 



m 

1.46 
0.24 
72 







-314-152i 


-290-150i 


-439-77i 


-383-12H 


a% 




-370-54i 


-340-80i 






a 5 




-55.7-31.4i 


-50-20i 


-45 


-27 


< 


-3? 


i0+(-46-106i) 


-380 


-575-77i 


-435-12H 




-3 


80+(-71-41i) 


-380 










1.25+0.3i 




0.5-1. 3i 


-0.89-2.73i 


a 8 


3. 


8+(-0.1-0.5i) 




4.6-0.4i 


3.3-0.9H 


a 9 




98.4+1.471 




-94-1.3i 


-93.9-2.7i 


aio 




39.3+7.23i 




-14-0.4i 


0.32-0.90i 
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For the neutral B decaying into CP eigenstate /, i.e., f — f, the effects of -B — B° mixing 
should be taken into account in studying CP asymmetry. Thus the CP asymmetry is time 
dependent, which is given by [pT5| 

A C p{t) = A d J P cos(Amt) - 2Im ^ CP h in{Amt), (32) 

where Am is the mass difference of the two mass eigenstates of neutral B mesons, and Aq P 
is the direct CP asymmetry defined in eq.(j3~lD with replacement of B + — > B° and B~~ — > B°, 
respectively. The parameter Xqp is given by 

Xcp - v tt v,i(f\H etf \B°y (33) 

With the above parameters and formulae, we get the branching ratios 



BK(B° - 


-> 7T + 7T ) 


= 7.55 x 10" 


■ 6 |e" 


-^ + .18e i8 -°T, 




- 7r°7r°) 


= 4.3 x 10~ 8 


|e" ; 


^ + 1.19e- a32 °| 2 , 


BK{B- - 


-> Tl\~) 


= 4.73 x 10" 


■ 6 |e" 


-^ + 0.05e- iai °| 2 , 


BTZ(B° d - 


■* K°7l°) 


= 4.06 x 10" 


■ 9 \e 


-^ + 31.9e i34 °| 2 , 


BK(B° d - 


■ K~Ti + ) 


= 5.12x10" 


"V 


-^ + 5.23e- 4l72 °| 2 , 


BR(B- - 


» K tt°) 


= 2.91 x 10" 


■ 7 \e 


-*T + 5.78e- fl68 T, 


BK(B~ - 




= 4.08 x 10" 


■ 9 \e 


-^ + 55.1e- ilL3 T 



(34) 



If we generally express eq.flSl]) as BIZ = A(e * 7 + ae tS ), then the direct CP asymmetry in 
eq.([H]) can be relevantly expressed as 

A dir = 2Qgm 7 f35) 

CP 1 + a 2 + 2acos5cos-f 



Using the above equation, the numerical results for the direct CP asymmetry are obtained 

5.0% 
1.03 + 0.36cos7* 



Af P {B -> tt+tT) = — sm 7 , 



A%(B -> ttV) = — sm 7 , 

CPK ' 2.42 - 1.59cos7 

Afp{B -> 7r°7r T ) = -1.7 x 10" 4 sm 7 , 

^(fi^A'V) = 3.5%szn 7 , (36) 
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A d j P (B ^ K^it*) = — — j— siwy, 



1.46 

28.4- 10.4cos7~ 



2 40 

cpy 1 34.4- 11.3cos7 

A d S r P {B K°tt t ) = -0.7%sm 7 . 

As is shown in Eq.(^), the strong phases are different by decay channels. We can also see 
from eq.(|36|) that the direct CP violation in B — > 7T -ir T is neglectably small. The direct CP 
violation in B — > 7r + 7r _ , tt /^, K°ir , K T ir° and -fr^Tr 1 * 1 are only at a few percentage level. The 
large CP violation effect may be expected in B — > 7r°7r° decays. However, it would be remained 
undetectable before the running of the next generation B factories, for example, LHCB, due to 
its very small branching ratios(~ 10 -7 ) and its two neutral final states. 

Recently, the CLEO collaboration has made first observation of the decay modes B — > 7i + ir~ , 
B —> K°ir° and B — > I^ir and also updated the decay modes B — > K ± ti t and B — > K°ti ± as 
follows H 



BK{B d - 


-> 7T + 7T ) 


= (4.3±^;1 ± 0.5) x 10" 6 ; 


BK(B U - 


-> 7T 7T ± ) 


< 12.7 xlO" 6 , 


BK(B d - 


-> ^V) 


= (14.6±I;?±I;|) x 10" 6 , 


BK(B d - 


► ^TT*) 


= (17.2±|;| ± 1.2) x 10" 


BTZ(B U - 


■* ^TT ) 


= (H.6±l;°±l;|) x lO" 6 , 


BTl(B u - 




= (18.2±t:S ± 1.6) x lO - 



(37) 



To compare with the data, we plot the CP averaged branching ratios for those modes as a 
function of 7 in Fig.2. Our results are plotted as curves and the CELO data are displayed as 
horizontal lines ( thicker lines for center value, thin lines represent error bars at 2a level). The 
horizontal line in Fig. 2. 7 is the upper limit of the decay mode. 

We find that the observed branching ratios of those decay modes can be well accom- 
modated within the QCD improved factorization approach of RefH except the decay mode 
B — > K°7T°. As shown in Eq . (|T9|) , the first term with F B ^ K and the second term with F B ~* 7T 
are dis constructive which reduces the amplitude of MiB — > K°7r°) much smaller than that 
of other B — > irK decays. As it is argued in Ref0, [T6j, in the present theoretical frame- 
work, the final state interactions are computable and identical to the imaginary part of the 
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Fig.2.1, Br(B -> K^ir^) vs 7 
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Fig.2.3, Br{B K^tt ) vs 7 
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Fig.2.5, Br(B -> vr^) vs 7 
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Fig.2.2, Br(B -> K°tt°) vs 7 
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Fig.2.4, Br{B -> K°tt^) 



vs 7 
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Fig.2.6, Br(B -> 7r°7r°) vs 7 
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Fig.2.7, Br(B -> vr^vr ) vs 7 

Figure 2: CP-averaged BR(B — > -kit, Kit) as a function of 7 are shown as curves for F B ^ 1T = 0.3 
and l^b/V^fel = 0.08 (in units of 10 -6 ). The BR measured by CLEO Collaboration are shown by 
horizontal solid lines. The thicker solid lines are its center values, thin lines are its error bars or the 
upper limit. 
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amplitude which is generated by the hard scattering amplitudes. In this paper, we find the 
strong phase appears not large enough to change the two sub- amplitudes of MiB — > K°7r°) 
to be constructive. Our results agree with that in Ref . fT3|, [17], [Tj| [19|] where the decay rate of 
B — > K°7T° is also estimated to be small. 

The CLEO observations have motivated many theoretical studies of those decay modes 
using different approaches [11], [12|, [17|, [Tj| |20fl . In Refs.[l8|, [H], it is suggested that 7 > 90° 
is required to interpret the CLEO data. However, the global CKM fit has given the constraint 

The comparison between our results and CLEO data|| implies 



7 < 90° at 99.6% C.L. 



120° < 7 < 240° which arises from the constraint by Br(B -^71 7r + ). The observed BriB 



7T 7T 



is smaller than many theoretical expectations. Negative cos'-f is needed to suppress 



the theoretical estimations as it is suggested in Ref.[TS[]. The decay rate of B — ► 7t~tt + can 
be also suppressed by using smaller form factor F B ^ n (0) and/or smaller \V u b/V c i,\. However, 
it would be very hard to account for the large decay rates of B — > Kit modes in this case. 
For those reasons, it might be difficult to solve the controversy between the global CKM fit 
and the mo del- dependent constraints from the charmless decays B — > Ktt, nix within the QCD 
improved factorization approach. 



4 Summary 

We have studied B A^tt ± B -> K°n° B -> K*ir° B -> K°n^ B -> 7r T 7r ± B -> ir°n° and 
B — > 7i T 7i° decays, in QCD improved factorization approach. 

The strong penguin contributions (Fig.l.e,f) are discussed in detail and found to be small 
because of the cancellations between them. The most important power corrections to these 
chiral enhanced terms(i.e., a^) are identified and found to be free of infrared divergence. With 
the choice of twist-3 DA (f) p {x) = 1, the gets a large imaginary part and its real part is 
enhanced by 10 ~ 20%. The other NLO coefficients aj also acquire complex phases from 
the hard scattering as depicted by Figl.(a) ~ (e) which are shown by the function g(x) and 
G(s, x) in Eq.(P^). We can see that g(x) is a new source of strong phase besides G(s, x) of the 
well known BSS mechanism^]. Compared to the Naive factorization, the strong phases are 
estimated reliably without the arbitrariness of gluon virtuality k 2 within the QCD improved 
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factorization formalism[|J. The strong phase due to the hard scattering in the decay modes are 
found to vary from 0° to 172°, depending on the decay mode. In the decays B — > 7r°7r°, K ± tt t 
and K ± 7i°, the strong phase are found to be as large as 100° < 5 < 180°. In other decay modes, 
the strong phases are rather small. 

The predicted branching ratios of B — > tiK and B — > -k^-k^ decay modes are in good 
agreement with the experimental measurement by the CLEO Collaboration except the decay 
B — > fT 7r°. The most serious constraint on the weak angle 7 comes from the small experimental 
value of Br(B — > 7r~7r + ) which implies 120° < 7 < 240°. We found that it is hard to solve the 
controversy between the constraints on 7 from the global CKM fit and the estimations of the 
charmless decays B — > Kit, tttt. The CP violation effects in B — > 7r°7r =F is neglectably small. 
The direct CP violation effects in B — > ir + 7i~ , tt°K t , K°ti , K^ti and K^ir^ 1 are only at a few 
percentage level. The large CP violation effect may be expected in B — > 7r°7r° decays. 

Note added: After finishing this work, we find Ref.|25[ also discussed B — > Kir and irn 
decays with a similar method, and Ref . P6| compared different approaches. 

Addendum After the paper was sent for publishing, the BARBAR Collaboration report 
their measurment of branching ratios for charmless B decays to charged pions and kions|27||: 
B(B° -> t^ttT) = (9.3±i;|±i;l) x 1(T 6 and B(B° -> K^) = (12.5±|:g±i;f) x 10~ 6 . Our 
predictions agree with the BARBAR data very well. We note that positive cos 7 is favored if 
the BARBAR data is taken as guide. 
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